Ab-initio Gorkov-Green's function calculations of open-shell nuclei 
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We present results from a new ab-initio method that uses the self-consistent Gorkov Green's 
function theory to address truly open-shell systems. The formalism has been recently worked out 
up to second order and is implemented here in nuclei for the first time on the basis of realistic 
nuclear forces. We find good convergence of the results with respect to the basis size in 44 Ca and 
74 Ni and discuss quantities of experimental interest including ground-state energies, pairing gaps 
and particle addition/removal spectroscopy. These results demonstrate that the Gorkov method is 
a valid alternative to multireference approaches for tackling degenerate or near degenerate quantum 
systems. In particular, it increases the number of mid-mass nuclei accessible in an ab-initio fashion 
from a few tens to a few hundreds. 

PACS numbers: 21.10.-k, 21.30.Fe. 21.60.De 
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Introduction. - The reach of ab-initio nuclear struc- 
ture calculations has extended tremendously over the last 
decade. Methods such as coupled-cluster (CC) [ill, in- 
medium similarity renormalization group (IMSRG) Q or 
Dyson self-consistent Green's function Q (Dyson-SCGF) 
have accessed medium-mass nuclei up to A ~ 60 on the 
basis of realistic two-nucleon (2N) and phenomenologi- 
cal three-nucleon (3N) 0] forces. In their current imple- 
mentations, such methods are however limited to dou- 
bly closed (sub-)shell nuclei and their immediate neigh- 
bors [5, 6]. As' one increases the nuclear mass, longer 
chains of truly open-shell nuclei connecting isolated dou- 
bly closed-shell ones emerge and cannot be accessed with 
existing approaches. Many-body techniques that could 
tackle genuine (at least) singly open-shell systems would 
immediately extend the reach of ab-initio studies from a 
few tens to several hundreds of mid-mass nuclei. It is the 
aim of the present Letter to propose one manageable way 
to fill this gap. 

Typically, open-shell systems can be dealt with via 
multireference schemes such as, e.g., multireference 
CC or configuration interaction techniques based on 
microscopic one- and two-body components of valence- 
space interactions 0, [f|. Alternatively, we presently 
wish to keep the simplicity of a single reference method. 
This requires however, in any of the approaches men- 
tioned above, to formulate the expansion scheme around 
a reference state that can tackle Cooper pair instabili- 
ties, e.g. to build the correlated state on top of a Bo- 
goliubov vacuum that already incorporates static pair- 
ing correlations. The simplicity of the single reference 



* vittorio.soma@physik.tu-darmstadt.de 
t C.Barbieri@surrey.ac.uk 

* |thomas.d uguet@c ea.fr| 



method can thus be kept at the price of breaking the 
symmetry associated with particle number conservation. 
The associated (small) contamination of the results that 
arises in finite systems eventually calls for the restora- 
tion of the broken symmetry Q . Notice that the use of 
Bogoliubov-based many-body methods to deal with near- 
degeneracies and non-dynamical correlations has recently 
been imported to quantum chemistry and proven to be 
extremely powerful [Ioj |. Further extension to calculate 
affinities and ionization energies would require an elec- 
tron attachment/removal formalism as the one we em- 
ploy. 

The present work discusses the first results obtained 
by extending Dyson-SCGF theory to the Bogoliubov al- 
gebra [ll|, i.e. it carries out the ab-initio Gorkov-SCGF 
formalism 12] in finite nuclei for the first time. A spe- 
cific benefit of such a method is to provide a way to 
understand microscopically and quantitatively the pro- 
cesses responsible for the superfluid character of atomic 
nuclei [13]. As a first step, normal and anomalous self- 
energies are truncated at second-order on the basis of 
2N interactions only. This constitutes a Kadanoff-Baym 
^-derivable approximation, i.e. the exact fulfillment of 
conservation laws is automatically ensured [l4j . As such, 
the method involves dressed propagators and is thus in- 
trinsically non perturbative. In the short-term future, 
the objectives are to incorporate 3N interactions into the 
frame and to generalize state-of-the-art Faddeev random- 
phase approximation (FRPA) truncation scheme [TH, [l6[ 
to the Gorkov context. 

Below, we present converged proof-of-principle calcu- 
lations of open-shell 44 Ca and 74 Ni nuclei. We employ 
a next-to-next-to-next-to-leading-order (N 3 LO) 2N chi- 
ral interaction [l7[ (A x =500 MeV) complemented by the 
Coulomb force. The resulting isospin-symmetry breaking 
interaction is softened using free-space similarity renor- 
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malization group (SRG) techniques [18j down to a mo- 
mentum scale of A = 2.0 fm _1 . After providing basics of 
Gorkov-SCGF formalism, we illustrate the binding en- 
ergy convergence with respect to the size of the har- 
monic oscillator basis used to expand the many-body 
problem. Then, various observables of experimental in- 
terest, i.e. ground-state binding energy, radii and pair- 
ing gaps, as well as adjacent isotopes' spectroscopy are 
discussed. Eventually, the effective neutron shell struc- 
ture [H GJ, Hi is displayed. 

Method. - Results displayed in the present Letter 
strictly rely on the formalism and the numerical imple- 
mentation detailed in Refs. JjTjJ and [2l|, respectively. 
Given the intrinsic Hamiltonian H- mt = T + V — Tcm, 
Gorkov-SCGF theory targets the ground state |^o) of 
the grand-canonical-like potential fi = Hj nt — ji A, where 
fi is the chemical potential and A the particle-number op- 
erator, having the number A = ( , J , o| J 4| , i , u } of particles on 
averag^J The complete one-body information contained 
in |\l/o) is embodied in a set of four Green's functional 
G" (uj) known as Gorkov propagators [l2| . Their matrix 
elements read in the Lehmann representation as 
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The poles of the propagators are given by ujk = fi/c ~ fioi 
where the index k refers to normalized eigenstates of fi 
that fulfil 



(2) 



The residue of G 99 ' (uj) associated with pole cot relates to 
the probability amplitude Uj. (Vfc) to reach state \^k) by 
adding (removing) a nucleon to (from) ^o) on a single- 
particle statcQ. 

Self-consistent, i.e. dressed, propagators are solutions 
of Gorkov's equation of motion 



T + S n (w) 
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1 Any consideration associated with A=N+Z applies in fact sepa- 
rately to the number of protons Z and to the number of neutrons 
N. 

2 Vectors and matrices defined on the one-body Hilbert space "Hi 
are denoted as bold quantities throughout the paper. 

3 The component of vector associated with a single-particle 
state a is denoted by Correspondingly, the component as- 
sociated with the time-reversed state a (up to a phase r) a ) is 
denoted by U% EH. 



whose output is the set of (U, V)fc and UJk- Equation ^ 
reads as an eigenvalue problem in which the normal 
[£ n (w) and £ 22 (w)] and anomalous [£ 12 (w) and S 21 (w)] 
irreducible self-energies act as energy- dependent poten- 
tials. Eventually, the total binding energy of the A-body 
system is computed via the Koltun-Galitskii sum rule [22[ 

E ° = hi J ct dLJ TTHi [g11(w) [t + (/i + w) 1]] • (4) 

Separation energies between the A-body ground state 
and eigenstates of A ± 1 systems are related to the poles 
UJk through 



E± ee ,! ± ujk = ± [{* k \H illt \y k ) - (* |^i„t|*o) 
TM[<*fc|^l*fc)-(A±l)] , 



(5) 



where the error associated with the difference between 
the average number of particles in state l^fe) and the 
targeted particle number A ± 1 is taken care of by the 
last term of Eq. (|5|). Spectroscopic factors associated 
with the direct addition and removal of a nucleon are 
defined as 



SF+ ee Tr ni 



U fc U 



and SF fc -=Tr^[V£V?].(6) 



In open-shell nuclei, the odd-even staggering of nuclear 
masses is a fingerprint of pairing correlations and of- 
fers, through finite odd-even mass difference formulae, 
the possibility to extract the pairing gap. The most 
commonly used [23| three-point-mass difference formula 

An (A) equates the pairing gap with the Fermi gap 
in the one-nucleon addition/removal spectra E^ , e.g. 

Al 3) (A) ee (-1) A [E+ - Eq]/2. One-body observable 
such as mass or charge radii can be easily computed from 
G n (uj) [nj]. Moreover, effective single-particle energies 
(ESPE) introduced by Baranger as centroids e™ nt of one- 
nucleon addition and removal spectra E^ can be natu- 
rally computed in the present context [llj . Last but not 
least, the normal self-energy S 11 (w) is identified with the 
microscopic nucleon- nucleus optical potential I24j. 1291 , al- 
lowing for the computation of scattering states [26| . 

Proceeding to an actual Gorkov-SCGF calculation re- 
quires to truncate the diagrammatic expansion of the four 
self-energies S" (uj). As opposed to perturbation the- 
ory, the expansion involves skeleton diagrams expressed 
in terms of dressed propagators solution of Eq. Such 
a key feature of self- consistent Green's function methods 
allows the re-summation of self-energy insertions to all or- 
ders and makes the method intrinsically non-perturbative 
and iterative. In the present application, self-consistency 
is limited to the static, i.e. energy-independent, part [27] 
H" (oo) of the self energy. This constitutes the so- 
called "scO" approximation that grasps the dominant 
fraction of self-consistency effects at a tractable numer- 
ical cost pjl [2l|. At first order in vacuum interactions, 
Eq. ([3]) reduces to an ab-initio Hartree-Fock-Bogoliubov 
(HFB) problem with energy-independent normal and 
anomalous self-energies accounting for Hartree-Fock and 
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FIG. 1. (Color online) Binding energy for 44 Ca (upper panel) 
and 74 Ni (lower panel) as a function of the harmonic oscillator 
spacing frw and for an increasing size N max = max (2n + I) 
of the single-particle model space. Results are from (scO) 
second-order Gorkov-SCGF calculations. The inserts show a 
zoom on the most converged results. 

Bogoliubov diagrams, respectively. In the present ap- 
plication, both first- and second-order diagrams are in- 
cluded Hj. 

Gorkov's equation (j3|) can be transformed into a 
energy-independent eigenvalue problem of larger dimen- 
sionality and solved iteratively. The algorithm in- 
cludes a readjustment of the chemical potential to ensure 
that proton and neutron numbers are correct in average. 
The most dramatic aspect of the implementation consists 
of handling the increased dimensionality, i.e. the growing 
number of poles in G" (ui) [Eq. (JTJ] with iterations, by 
means of a Lanczos algorithm. All numerical aspects will 
be reported on in a forthcoming publication [21] . 

Results. - Figure Q] displays the binding energy of 44 Ca 
and 74 Ni as a function of the harmonic oscillator spac- 
ing hu and for an increasing size N max = max (2n + I) 
of the single-particle model space. The convergence pat- 
tern obtained here on the basis of a soft 2N interaction 
is similar to those generated for doubly-closed shell nu- 
clei with currently available ab-initio methods. Overall, 
convergence is well attained for N max = 13. In 44 Ca, 
going from N max = 11 to N max = 13 lowers the minima 
by just a few keV. Also, the binding energy calculated 
for N max = 13 varies by less than 200 keV over a wide 
range of Huj values. In 74 Ni, going from N max = 11 to 
Nmax = 13 yields an additional 600 keV, while scanning 
a large range of oscillator frequencies only changes the 
binding energy by about 1 MeV. 

Table U lists the results obtained for various observ- 
ables of interest in the ground state of 44 Ca and 74 Ni. 



The values quoted are extrapolated to infinite oscillator 
basis size using the method proposed in Ref. [28[ . At this 
point, results are mostly illustrative because of the lack of 
3N forces. The latter play a key role in the saturation of 
nuclear matter such that omitting it generates too much 
binding and too small nuclei when using soft 2N inter- 
actions [l8|]. The neglect of 3N forces also induces too 
small pairing gaps as a result of the too low density of 
states in the nucleon addition and removal spectra (see 
below). It is our short-term objective to add 3N forces 
to the present theoretical scheme. 
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FIG. 2. (Color online) One-neutron addition and removal 
spectral strength distributions in 44 Ca obtained from second- 
order (scO) Gorkov-SCGF calculations. For each final state 
in 43 Ca (left to the dashed line) and in 45 Ca (right to the 
dashed line), the spectroscopic factor is plotted as a function 
of its separation energy to the ground state of 44 Ca. Energies 
above MeV correspond to n+ 44 Ca scattering states [261 ]. 
Final states with different J 71 values are separated for clarity. 
Results correspond to the minimum of the convergence plots 
shown in Fig. [T] Although center of mass motion is subtracted 
by using Hi nt , the variation of that correction going from A 
to Ail is neglected. The associated error is small in such 
medium- mass nuclei [2111 . 



Figure [5] displays one-neutron addition and removal 
spectral strength distributions (SSD) in 44 Ca. Results 
are shown over a large range of final states in 43 Ca and 
45 Ca characterized by spectroscopic factors as small as 
2.10 -3 (i.e. 0.2%). One observes a fragmentation of the 
spectroscopic strength that is characteristic of correlated 
many-body systems. Overall the pattern is similar to the 
one found in doubly- magic nuclei [3j . Close to the Fermi 
energy, however, one notices a feature that is unique to 
open-shell nuclei, i.e. the 7/2~ strength is equally frag- 
mented into additional and removal channels, which re- 
sults in the fact that both 43 Ca and 45 Ca ground-states 
have angular momentum and parity J 71 ' = 7/2~. Such 
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a fragmentation reflects static pairing correlations that 
manifest themselves as a result of emerging degenera- 
cies in the ground state of open-shell nuclei. It is the 
main strength of Gorkov-SCGF theory to explicitly han- 
dle such degeneracies and resulting pairing correlations. 

The right column in the upper panel of Fig. [3] supplies 
a zoom of Fig. [5] around the Fermi energy for states with 
spectroscopic factors larger than 10 _1 (i.e. 10%). The 
left column provides the same quantities for first-order 
(i.e. HFB) calculations. Last but not least, the center 
column displays effective single-neutron energies. The 
same information is provided for 74 Ni in the lower panel 
of Fig. El 

The main fragmentation of the strength is absent from 
first-order calculations, i.e. it is due to dynamical corre- 
lations that come in at second order and that are qualita- 
tively the same as for closed-shell nuclei. Contrarily, the 
fragmentation of the strength in the vicinity of the Fermi 
energy into two peaks of (essentially) equal strength is 
qualitatively accounted for at first order and thus relates 
predominantly to static pairing correlations. Quantita- 
tively, the energy spacing between the two low-lying 7/2~ 
states in the SSD of 44 Ca is increased by second-order ef- 
fects. This is in contrast to the behavior of 74 Ni where 
the separation of the low- lying 9/2+ states is instead de- 
creased. Given that such a spacing is equal to (twice) the 
pairing gap, one concludes that the coupling of Cooper 
pairs to non-collective fluctuations may already affect 
pairing correlations in open-shell nuclei in either way. A 
detailed study of such a feature is left to a forthcoming 
publication^. Back to the overall spectrum, one observes 
that the position of the dominant peak of a given J n value 
is significantly modified by second-order effects such that 
the corresponding spectrum is more compressed than at 
first order. However, it is still significantly too spread 
out compared to experiment due to missing 3N forces 
and the lack of coupling to collective fluctuations. 

Effective single-particle energies recollect the frag- 
mented strength [ill El S3 from both one-nucleon ad- 
dition and removal channels. Many-body correlations 
are largely screened out from ESPEs, which picture the 
averaged single-nucleon dynamics inside the correlated 
system. Two different features are identifiable in the 
ESPE spectrum ejf nt when compared to observable one- 
nucleon addition and removal spectra Est- The ESPE 
e i/"/2 ( e igg/ 2 ) l° ca ted at the Fermi energy recollects 
the strength of the two equally important 7/2~ (9/2 + ) 
states. Other ESPEs recollect the strength of a low- 
lying dominant peak and of a highly fragmented strength 
distributed at higher excitation energies such that they 
move away from the Fermi energy to closely match first- 



4 A quantitative treatment of nuclear superfluidity through ab- 
initio approaches requires to treat the coupling of the Cooper 
pair to collective density, spin and isospin fluctuations [l3, 29]. 
In the present context, this necessitates the implementation of 
the (generalized) FRPA expansion scheme |15|. |16( . 
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FIG. 3. (Color online) Left: One-neutron addition and re- 
moval spectral strength distribution obtained from first-order 
(HFB) Gorkov-SCGF calculations. Right: same as left panel 
for second-order (scO) calculations. Center: Baranger ESPEs 
reconstructed from second-order (scO) Gorkov-SCGF calcula- 
tions. Upper panel: 44 Ca. Lower panel: 74 Ni. 

Conclusions. - Wc have presented the results of the 
first-ever ab-initio calculations of medium-mass (truly) 
open-shell nuclei. Such calculations are based on the 
implementation of self-consistent Gorkov Green's func- 
tion theory on the basis of realistic nuclear interactions. 
Taking 44 Ca and 74 Ni as test cases, we have demon- 
strated the good convergence of the results with re- 
spect to the basis size and discussed several quantities 
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of experimental interests including ground-state energies, 
pairing gaps and particle addition/removal spectroscopy. 
Such calculations increase the reach of ab-initio calcu- 
lations in the mid-mass region tremendously and are 
now being performed systematically over long isotopic 
and isotonic chains [21] . The short-term objectives are 
to incorporate three-nucleon interactions into the frame- 
work and to extend state-of-the-art Faddeev-random- 
phase-approximation truncation scheme from doubly- 
closed shell nuclei to open-shell nuclei, i.e to the present 



Gorkov context. 
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